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We explicitly evaluate the infinite series of integrals that appear in the "Anderson- Yuval" reformu- 
lation of the anisotropic Kondo problem in terms of a one-dimensional Coulomb gas. We do this 
by developing a general approach relating the anisotropic Kondo problem of arbitrary spin with 
the boundary sine-Gordon model, which describes impurity tunneling in a Luttinger liquid and in 
the fractional quantum Hall effect. The Kondo solution then follows from the exact perturbative 
solution of the latter model in terms of Jack polynomials. 
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It was shown some 25 years ago that the partition func- 
tion of the anisotropic Kondo problem could be expressed 
as a power series in the impurity coupling Even 
though the coefficients were given only as formal multi- 
dimensional integrals, this led to progress in the under- 
standing of the Kondo problem. For example, this en- 
abled renormalization-group calculations to be done and 
the phase structure determined ||. In subsequent years, 
various additional aspects of the problem were further 
understood by methods like the numerical renormaliza- 
tion group 0, expansions around the low-temperature 
fixed point || and by the Bethe ansatz Despite 
all this progress, the problem is still far from completely 
solved. For example, the finite-temperature resistivity 
(the aspect which motivated Kondo's analysis 0) has 
yet to be calculated exactly. Moreover, the coefficients of 
the perturbative expansion still could not be evaluated 
explicitly. 

This latter aspect highlights a general puzzle of the 
Bethe ansatz approach. While its results are non- 
perturbative, only by crude numerical solution does it 
provide much information on the small- coupling pertur- 
bative expansion (except when T — 0). The reason is 
that for T ^ 0, the Bethe ansatz gives a solution in 
terms of non- linear integral equations, and no technique 
is known to solve them systematically around the non- 
interacting fixed point. In this paper, we describe a new 
formalism which avoids these problems. We show how to 
"do" the Anderson- Yuval integrals, by expressing them 
in terms of infinite sums of ratios of gamma functions 0] 
which can be easily evaluated numerically. 

We derive a simple functional relation between the par- 
tition functions of the Kondo problem and another well- 
known one-dimensional quantum problem, the bound- 
ary sine-Gordon (BSG) model. The derivation utilizes 
some results of [jl(J , together with the realization that the 
BSG model can be considered as a particular anisotropic 
Kondo problem where the boundary spin is in a "cyclic" 
representation of the quantum group SU(2) q . The per- 
turbative coefficients of both of these partition func- 
tions can be expressed as the partition function of a 
one-dimensional classical "log-sine" gas with positive and 



negative charges which have logarithmic interactions. In 
the Kondo problem the charges must alternate in sign in 
space, while in the BSG model they may occur in any 
order. With unrestricted ordering, the partition func- 
tion can be evaluated J9[ by utilizing various properties 
of Jack polynomials [ |ll| . Since the functional relation 
gives the appropriate Kondo coefficients simply in terms 
of the BSG ones, it therefore allows their determination 
as well. 

We first review the BSG and Kondo perturbative ex- 
pansions. We then derive the (non-perturbative) relation 
between these partition functions, which gives one set 
of perturbative coefficients in terms of the other. Both 
sets have applications in a variety of physical problems. 
The Kondo problem of course is realized in various impu- 
rity compounds, and is also equivalent (in the anisotropic 
case) to the dissipative quantum mechanics of a particle 
in a double well |l2| . The BSG model describes tunnel- 
ing in a Luttinger liquid (with application to fractional 
quantum Hall edges) |l3||l4|], while in dissipative quan- 
tum mechanics it corresponds to an infinite number of 
wells g5|. 

The boundary sine-Gordon model describes a free bo- 
son <f>(<7, t) on the half-line a > with an interaction at 
the boundary a = 0. We study the problem at non-zero 
temperature T; in the path integral this corresponds to 
a system where Euclidean time (denoted by r) is a circle 
of periodicity 1/T. The bulk action is 
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while the boundary action is 



dr \{d T cj>) 2 + {d a <pf 



2A 



1/T 



dr cos [</>(er = 0, r)] 



(1) 



(2) 



The coupling A is not scale-invariant and interpolates be- 
tween Neumann (A = 0) and Dirichlet (A — > oo) bound- 
ary conditions on the boson. 

Defining as usual the partition function via the path 
integral as Z = f[d(f>] e s+Sn we introduce Z = 
Z(X)/Z(X = 0). We can rewrite this as a series by ex- 
panding out exp^Ss) in powers of A. The coefficients of 
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this series are then integrals of correlators in the free- 
boson theory. Using the bosonic propagator on the edge 
of a half-cylinder with Neumann boundary conditions, 



MO,TMoy)>=-20ln — sinTrT(T-r') 
(k is the inverse cutoff), one finds for example 



(e 



#(0,t)_-#(0,t')\ _ 



nT 



sin7rT(r — t ) 



-2fl 



Using Wick's theorem one finds the general correlators 
as ratios of various C(t,t') = |2sin7rT(T — r')| 2s . The 
partition function is the series 



Z = J2 x 2n Z 2n 



n=0 



(3) 
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where all indices run from 1 to 2n, and all integrals run 
from to l/T independently. This expression requires 
regularization for g > 1/2, but is well-defined otherwise. 
The sum (^) is the grand canonical partition function 
for a classical two-dimensional Coulomb gas of charged 
particles restricted to lie on a (one-dimensional) circle 
of circumference l/T. The charges are ±*/2g, and the 
whole gas is electrically neutral; x is the fugacity. 

Explicit series expressions were found for the Z 2n || by 
expanding the denominator of (Q) in terms of Jack poly- 
nomials |ll[] , and using orthogonality relations of these 
polynomials. Jack polynomials are indexed by a par- 
tition of an integer m = (mi,m2, "^((m))); where 

n%i > m 2 . • • > TO;( m )- Any m can be pictorially repre- 
sented by a Young tableau with l(m) rows and rrij blocks 
in the j^ 1 row. One finds 0] 
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where the sum is over all partitions m with l(m) < n. 
The sum converges only for g < 1/2, which is where the 
integrals in (^) are well-defined. 

For n = 1, the sum can be done explicitly (or the 
integral can be done without Jack functions) giving Z 2 = 
T(l — 2g)/r 2 (l — g). Notice that the divergence at g = 
1/2 (the well-known "free-fermion point") is just a simple 
pole, so we can analytically continue around it to study 
behavior for g > 1/2. In fact, it turns out that at the free 
fermion point, only the first coefficient in the expansion 
of the free energy (= —Tin Z) diverges making it possible 



to analytically continue all the Z 2n to g > 1/2 Q. The 
pole in Zi is a signal that there are logarithmic terms in 
the perturbative expansion at g = 1/2, which can in fact 
be computed exactly [ ^6[ . Thus this free-fermion point 
(which corresponds to the Toulouse limit of the Kondo 
problem) is pathological in some respects. 

The perturbative expansion for the Kondo problem 
found in [QJ^] is closely related. The Kondo problem is 
a three-dimensional non-relativistic problem, with free 
electrons antifcrromagnetically coupled to a single fixed 
impurity. By looking at s-waves only, we restrict to the 
radial coordinate and this becomes a one-dimensional 
quantum problem where massless fermions tpi (cr, t) (i — 
1, 2 are the spin indices) move on the half-line a > with 
a quantum-mechanical spin S a at a — 0. The interaction 
parameters are I Z ,I+ = I- '■> in the isotropic case I z = I+- 
The impurity action is 



l/T 



dT^t(0,T)S o o§V;(0,T) (6) 



where the a a are the Pauli matrices. By a well-known 
bosonization procedure j^j , this can be rewritten in terms 
of a free boson <f>(x,t) with bulk action (|l|) and 



Sb =h 



l/T 



dr (s*+e #( °' r) + &.c) , 



where we have absorbed I z in a redefinition of the cou- 
pling g. The anisotropy is thus parametrized by g, with 
g = 1 the isotropic case and g = 1/2 the Toulouse (free- 
fermion) limit. 

In the perturbative expansion, we get correlators of 
exp(±i</>) like before. The crucial difference arises from 
the S±. When the impurity has spin 1/2, we have 

S\ = Si = 0, and only terms of the form S+SS+S 

survive in the perturbative expansion. In the one- 
dimensional gas, this is the requirement that charges 
alternate in sign. We define the partition function of 
the spin- 1/2 Kondo problem as Zk — 2Zk(I+)/Zk{0), 
where the factor of 2 ensures that the entropy is In 2 in 
the non-interacting limit I + = 0. Thus 



Z K (X K ) =2 + Y J (xK? n Q2r, 
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The Q2n are the integrals in times 2(n!) 2 , but where 
the region of integration is restricted to be < t\ < 
t[ < t-x < . . . r' n < l/T. The periodicity of the integrand 
means that Z 2 — Q 2 , but the others are different. 

Our central result relates the partition functions of the 
Kondo and BSG models. Defining q = exp(i7T(7), we find 
that 

Z K [(q - q-^Zix) = Z(qx) + Z(q- X x). (7) 

One simple check is that the equation is consistent with 
the fact that Z(x) = 1 in the isotropic case q = — 1. 
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Plugging in the perturbative expansions gives the Kondo 
coefficients Q 2n in terms of the BSG coefficients Z 2n . For 
example 



4sin 2 (7rg)(34 = -4 cos 2 (7rg)Z 4 



(Z 2 f 



16sin 6 (7T3)Q 6 = sm 2 {?,itg)Z & + (Z 2 ) 3 sm 2 (irg) 
- (sin 2 (27Tff) + sin 2 (7r.g)) Z 2 Z 4 



(8) 
(9) 



Using the series (g), we can easily find values for the 
Z2n and hence the Q 2n for g < 1/2 by truncating the 
series and summing it numerically. For example, we have 
Q2 = ^2=r(l-2.g)/(r(l- 3 )) 2 ,and 



Values of CoefRcents 




g = 2/5 


g = 1/3 


.9 = 1/4 


z 4 


1.910750624 


0.837804224 


0.4644013099 


z 6 


1.088518710 


0.276783311 


0.0968299150 


Z 8 


0.439166887 


0.061847648 


0.0129159832 


Zio 


0.135465650 


0.010210054 


0.0012208002 


Qi 


0.982706435 


0.432237451 


0.2322006549 


Q 6 


0.291860092 


0.074496009 


0.0242074788 


Q s 


0.061852434 


0.008729195 


0.0016144979 


Qw 


0.010067801 


0.000757798 


0.0000763000 



As a check, some of the integrals were explicitly evalu- 
ated using Monte Carlo methods; the agreement is good. 
As discussed in Q , for g rational there are relations be- 
tween the numbers beyond those in (Q); for example for 
g = 1/4, Q 2n = Z 2 „/2 n - 1 . The knowledge of the first 
few coefficients Q 2n provides a very good numerical so- 
lution of the problem. The series can be extrapolated 
using Pade approximants and for instance the flow of the 
boundary entropy is reproduced fairly accurately . 

The proof of (Q) requires first of all the introduction 
of the "quantum-group" algebra SU(2) q , which is a one- 
parameter deformation of the SU(2) algebra |L7||. There 
are three generators S+ , S- and S z , with commutation 
relations 



qS^ 



1 



-s. 



[S+,S-] 
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q-q 



Like SU(2), the quantum-group algebra has represen- 
tations of any spin. A careful analysis shows that the 
general anisotropic Kondo problem of arbitrary spin is 
integrable only if the spins in (||) obey the SU(2) q alge- 
bra. For the isotropic problem q = — 1, the distinction 
is irrelevant because 5C/(2)_i is identical to SU(2). The 
distinction is also irrelevant for any q for spin 1/2 or 1, be- 
cause the the two-dimensional representation of SU(2) q 
is given by the Pauli matrices as for SU (2), and the spin- 
1 representation is also the same up to a rescaling of S+ 
and S-. Other representations differ, however, so one 
must be careful when discussing the anisotropic Kondo 
problem at higher spin. 



Following [ JL0[ , let us now introduce the "quantum 
monodromy operator" Lj(x), which is defined as 



Lj (x) = T exp 



cl/T 



Jo 
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Here the T indicates time-ordering, and the SU (2) q gen- 
erators are in the spin-j representation. The field 4>l is 
the left-moving component of the field </>; with Neumann 
boundary conditions we have 2^(0, t) = 0(0, r). The 
"quantum transfer matrix" Tj = Tr Lj is a fundamen- 
tal part of the Kondo problem: (Tj(x)) is identical to 
the spin-j anisotropic Kondo perturbative expansion. In 
particular, 

Z K (x) = (T 1/2 (x)). 

This follows from expanding (]l^) in powers of x (the 
factors q^ Sz cancel), and using the fact that the vacuum 
is an eigenstatc of Tj . 

As shown in jlo| , the Lj satisfy the Yang-Baxter equa- 
tion. This results in a number of remarkable properties. 
In particular, one finds that all of the Tj can be generated 
from spin- 1/2 via the relation jl0| 

Tj(q l ' 2 x)Tj(q- l ' 2 x) = 1 + Tj_ i(x)T j+ i (x) (11) 

or via jl(| 

TiW+^x^x) = T j+i {q*x) +Tj_ h (q-ix) (12) 

This gives enough information to relate the two partition 
functions Z and Zk (and do much more). When one 
expands Tj{x) for higher j in powers of x, one obtains 
integrals with all sorts of charge orderings, with weights 
depending on q (because 5+ and S- depend on q for 
representations other than spin 1/2, and because of the 
monodromy of the chiral fields). For example, S± = 

for spin 1, so + H h— appears at order x 6 in Ti, 

but + + H does not. Therefore, the unordered 

integrals Z 2n in (fj) can be constructed by summing over 
appropriate combinations of terms from the Tj{x). Since 
the relations ( |Tl| ) or ( fl2"| ) give all of the higher Tj{x) in 
terms of T1/2, this means that any Z 2n can in principle be 
expressed in terms of the Q 2n . The perturbative results 
(||) and (||) follow directly from ( . 

The quickest way to derive (u) is to introduce cyclic 
representations of SU(2) q . These occur when q is any 
root of unity, q l = ±1, and have no analog in ordinary 
SU{2). They are indexed by an arbitrary complex pa- 
rameter 8 and have t states |m), m = 0, . . . ,t — 1 with 
generators acting as 



S±\m) 



-<5±m 



-|m±l) S z \m) = m|m) 



3 



where states |m) and \m+t) are identified. Notice that in 
such a representation one can consider arbitrary powers 
of the generators. Referring to cyclic representations as 
spin-<5, one finds that all previous results for spin-j apply. 
In particular, relation holds with j — > S. In the 
formal limit where 5 — — iA with A >> 1 we have 



S± 



In this limit the commutator [S+, S-] becomes negligible, 
so the traces of all monomials with n generators S+ and 
n generators S- become identical, e.g. 



Tr(5+5_ 



t 



,1/2 TrgA 



(One might have been tempted to study lim^oo Tj where 
arbitrary powers of the generators also occur, but this is 
not Z because the traces do not behave in the right way.) 
Hence one way of expressing Z is 



Z(x) 



1 



■TvgA 



A»l 



(13) 



We substitute j — > 8 — —iA and x — * (q — q~ l )q~ s ~^ 2 x 
into (p"2|). Letting A >> 1 proves (0) for q any root of 
unity (rational g). The result follows for any \q\ = 1 
(real g) by continuity. The renormalization of the Kondo 
coupling in ([?]) can be checked through the identity Z 2 = 
Q2- This renormalization ensures that (Q) makes sense 
in the isotropic limit q — * — 1, where the Kondo partition 
function is a non-trivial object but Z(x) = 1. 

Having obtained the fundamental relation (^), higher- 
spin Kondo partition functions follow from fusion, using 
the relations (11) or ([L2|) . For instance one has 



1 ' 2 x) = 



(Tj [{q-q- l )x\)Z{q l / 2 x)Z{q 
Z{q 3 / 2 x)Z(q^ 2 x) + h.c. + Z(qV 2 x)Z{q- 3 / 2 x) 



(14) 



Other such relations give integrals of the form (Q) with 
any charge ordering. Slightly more complicated ordered 
integrals arise in many places, and we hope our results 
can be extended. One possibility is in the Keldysh calcu- 
lation of the non-zero-frequency noise in the BSG model 
|ji"8|| . Another is the quantity P(t,T) in the double- 
well problem of dissipative quantum mechanics, which 
is the probability that a particle at temperature T is in 
one well at time t given that it is localized in that well 
at t = @. We have found the curious result that 
P(iT, T) = Ti/2 — 1/2, but its significance is not clear. 

We have therefore succeeded in computing explicitly 
the Coulomb-gas integrals for the anisotropic Kondo 
problem. Our approach suggests a very deep struc- 
ture, where the integral equations of the Bethe ansatz 
can be solved in terms of Jack symmetric functions and 
quantum-group combinatorics. We hope that further 
studies will uncover more features of these relations. 
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